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Summary. These introductory lectures present a broad overview of the physics of 
f~ i . high parton densities in QCD and its application to our understanding of the early 

time dynamics in heavy ion collisions. 

■ Introduction 

Quantum Chromodynamics (QCD) is widely accepted as the fundamental 
^ I theory describing the behavior of hadrons. In QCD, hadrons are composed 

of elementary quarks and gluons which are often together labeled partons. 
Quarks and gluons are never directly measured due to the confining property 
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■ of the theory. However, their distributions inside a hadron can be probed 
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precisely in Deep Inelastic Scattering (DIS) experiments. It is seen from the 
HI and ZEUS experiments at HERA that the structure functions of the gluons 
QQ I and the sea quarks, which to leading order in QCD are related to their number 

■ densities, grow very rapidly [1] with decreasing values of a Lorentz invariant 

^ I kinematic variable xbj introduced by Bjorken. Again, to lowest order, this 

variable corresponds to the momentum fraction x of the hadron's momentum 
carried by a parton. Small x physics is the regime of high energies in QCD 
^ . and the physics of this regime exhibits many novel features that are not fully 

understood. 

Small X physics is interesting for a variety of reasons. Even if the momen- 
tum transfers squared are large enough such that, from asymptotic free- 
dom, the QCD coupling constant as is small, the explosive growth in the num- 
ber of partons with increasing energy makes the physics non-perturbative. The 
strongest electric and magnetic fields found in nature occur in this situa- 
tion. The answer to a fundamental question in QCD - can we calculate the 
total number of particles produced in a strong interaction at asymptotically 
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high energies - therefore requires that we understand the processes of parti- 
cle creation in the presence of such strong color fields. Particle production in 
this context is important for understanding a variety of striking but little un- 
derstood phenomena at high energies. These include, for instance, i)limiting 
fragmentation [2], where the rapidity distributions of the produced hadrons 
turns out to be independent of energy around the fragmentation region but 
possess non-trivial features in the central rapidity region, ii) the unusually 
large fraction of diffractivc final states in DIS where no particles arc produced 
in angular regions relative to the scattering plane called "rapidity gaps" [3] 
and iii) the phenomenon of "geometrical scaling" , where cross-sections appear 
to scale as a dimcnsionlcss function of the momentum transfer squared of the 
probe relative to a dimensionful dynamical scale in the hadron [4]. Small x 
physics is also of crucial importance in understanding the formation of bulk 
QCD matter in high energy heavy ion collisions and its possible thermalization 
to form a Quark Gluon Plasma (QGP). Finally, the hope is that understand- 
ing the properties of strong fields at small x will provide some insight into 
confinement in QCD and its role in high energy scattering. 

These questions can be addressed in a weak coupling effective field theory 
formalism [5, 6] which describes the properties of high energy wavefunctions 
as a Color Glass Condensate (CGC). In this approach, the degrees of freedom 
in the high energy nuclear wavefunction are divided into static light cone color 
sources at large x and dynamical gauge fields at small x which are coupled to 
these static color sources. Because the scale between the two sorts of degrees 
of freedom is arbitrary, and because physics cannot depend on this scale, 
one obtains a renormalisation group description in rapidity which arises from 
successively integrating out dynamical fields at one scale and absorbing them 
into color sources at the next. 

Efforts to understand the rich non-perturbative phenomena of high energy 
QCD are the subject of these lectures. The first lectmc begins with DIS and 
describes high energy scattering in DIS in the well understood Bjorken limit 
and the much less understood Regge-Gribov limit of QCD. We show that par- 
ton saturation arises naturally in the latter limit. The effective field theory 
formalism of CGC, which incorporates the physics of parton saturation, is 
described next. We then discuss color dipole models that incorporate simply 
both the non-linear dynamics of saturation and the linear dynamics of per- 
turbative QCD in DIS and hadron scattering. In the second lecture, we apply 
the CGC approach to treat ab initio heavy ion collisions at high energies. 
When two sheets of colored glass collide, they form bulk QCD matter called 
the Glasma in heavy ion collisions. We describe the properties of the Glasma, 
probes of its dynamics and the possible fast thermalization of the Glasma into 
a Quark Gluon Plasma. 
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1 Deep Inelastic Scattering 

Deep Inelastic Scattering (DIS) involves the scattering of a high energy lepton 
off a hadronic target in which the energy and momentum transfer of the lepton 
is measured experimentally. DIS is essentially a two step process where the 
lepton emits a virtual photon in the first step which then interacts with the 
hadron in the second step. Depending on the energy transferred, the proton 
can break up and new particles are created. This provides a clean environment 
to study the structure of hadrons at high energies since one of the interaction 
vertices in the two step process is completely known. The DIS process is 
schematically illustrated in fig. f . The cross section of the process is 



dxdQ'^ xQ* 




Fig. 1. Left. The general setup of DIS. Right. Inside the hadron via the Parton 
model. In the impulse approximation the partons inside are non-interacting on the 
interaction time-scale of the virtual photon probe. 



Light Cone coordinates : It is useful at this point to intro- 
duce light cone (LC) coordinates which are very useful to dis- 
cuss high energy scattering. Let z denote the longitudinal axis of 
collision. For an arbitrary 4-vector = {a* , a'^ , , a^) , the LC 
coordinates are defined as: 

j_ + _ a* — a^- , x v\ 

= —-j=-,a = — ^=— , aj_ = (a , a^). (2) 

In particular a;+ = (t + z) / is the LC "time" and ^ {t — 
z)/\/2 is known as the LC "longitudinal coordinate". 
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Here e and h denote the initial electron and hadron state and x the final 
hadronic state, a^m is the usual QED coupling constant. The four vectors of 
the electron and proton are shown in fig. 1. If denotes the (space-like) 4- 
momenta of the exchanged photon, then the energy transferred to the hadron 
is 

= -q^q^ = -{K- k'^f = ^E,E'^ sin^ (|) . (3) 

The Q'^ is a measure of the resolution power of the probe. and E'^ are the 
initial and the final energies of the lepton and 9[, is the lepton scattering angle 
in the center-of-mass(COM) frame. The variable y provides a measure of the 
inelasticity of the collision and is defined in the following frame invariant way 

P-q 1 -^e 2^e 

y = TTT = 1 - cos^ — (4) 
^ P.k Ee 2 ^ ' 

The X here is the Bjorken variable Xb^ whose frame invariant definition is 

^B, = ^ (5) 

At very high energies, wc obtain .tbj ^ Q'/s for a fixed y. This clarifies why, 
for a fixed Q^, small x physics is the physics of high energies in QCD. As we 
shall see shortly, Xbj is related to the momentum fraction of the struck quark. 

F2 is the structure function that, at leading order in QCD, gives the quark 
+ antiquark distributions in a proton, while the longitudinal structure func- 
tion Fl is a measure of the gluon momentum distribution. These quantities 
can be independently extracted from the DIS cross section by varying x, 
and the center of mass energy of the collision. The most striking feature ob- 
served in early DIS experiments was the flatness of the structure function F2 
over a wide range of at which the experiments were first performed (see 
figure 2). The apparent scale invariance of the structure function gave the 
needed experimental support for the hypothesis of point-like, weakly inter- 
acting partons inside hadrons. The dynamical model that was formulated by 
Feynman [7] to understand this phenomena is called the parton model. 



1.1 The Bjorken limit, the Parton model and pQCD 

The Bjorken limit, in which the results of the DIS experiments were correctly 

interpreted, is attained when the center of mass energy s 00 and the energy 
transferred in the collision ^ 00 keeping the ratio xbj = Q^/s = fixed. At 
these very high energies, the rest mass energy of the hadronic target M can be 
neglected compared to its longitudinal momentum P+. In the parton model, 
formulated in the "infinite-momentum" frame, the hadron can be thought 
of as a collection of "quasi-free" partons which are nearly on-shell excitations 
carrying some fraction xp of the total hadron momentum. Moreover the entire 
momentum of the hadron can be thought to be longitudinal. This picture is 



Two introductory lectures on high energy QCD and heavy ion coUisions 



5 



.J< 10 



do' 



10 



10' 



10' 
10^ 



10 



-0,000063 

x-0.000102 
x-0, 000162 
^=0.000253 
x=0.0004 
x=0,0005 
\=0, 000532 
. x=0.0008 



Prolo-i 

* HI 

* ZEUS 

* BCDMS 
O E665 

□ NMC 

A SLAG 



x=0, 00102 
=0.0013 
-0.00161 
x.0.0021 
X. 0.00263 
x-0.0032 
„. x.0.005 



-0.008 

x-0.013 
t x-0.021 
, * x=0.032 
. . . x-0.05 




-0.18 
x-0.25 

x=0.4 



=0.85 (i,= 1) 



10 1 10 10 10 10 10 10 

Q^(GeV^ 

Fig. 2. Bjorken Scahng: Independence of structure functions vs in proton. Cour- 
tesy PDG. 



essentially true since the interactions between the partons are highly time- 
dilated. In the "impulse approximation" when the virtual photon strikes a 
parton, the other partons form the spectators without interacting with the 
struck quark or among themselves. (See fig. 1). This demands 

{xfP + qf = ~ 0; xp ^ = Xbj , (6) 

which confirms our interpretation of Xbj ■ 

The naive parton model also predicts Fl to be zero. This result is called 
the Callan-Gross relation and it provides strong evidence that the partons 

probed by the virtual photon arc spin- 1/2 objects. It was also realized that the 
hadron might also contain an infinite sea of light qq pairs, called appropriately 
sea quarks, as opposed to the valence quarks that carry the net baryon number 
of a hadron. However it was experimentally found that, at a scale of 0(1 GcV), 
the proposed valence and sea quark distributions could only account for about 
50% of the total momentum in a proton. This necessitated the existence of 
other partons (gluons) which in turn explained the puzzle that Fl was not zero 
experimentally. With further experiments, it was also demonstrated that the 
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Bjorken scaling was only approximately true and there were scaling violations 
at lowcr-x. This did not have an explanation within the naive parton model 
at all. 

These experimental results set the stage for the formulation of the theory 
of strong interactions in terms of QCD which had quarks and gluons. Most 
important though, in tying all the pieces together, was the (Nobel) discovery 
of Gross, Wilczek and Politzer that the theory had a /3 function with negative 
sign- thereby indicating that the coupling asymptotically goes to zero at high 
Q^. Perturbative QCD (pQCD) therefore showed that in the Bjorken limit the 
naive parton model is indeed a good approximation. It went further to give 
a quantitative explanation of the log scaling violations that were observed 
in the experiments. The key idea to the solution of this puzzle lay in the 
assumption of neglecting the transverse momentum of a quark. In fact, a 
quark can emit a gluon and acquire a large transverse momentum pT, on time 
scales comparable with that of the probe, with a probability proportional to 
asdp^/p^ at large pt where as is the strong coupling constant. On integrating 
this till the kinematic limit of Px ^ Q^, contributions proportional to InQ^ 
are obtained which are precisely the scaling violations. 

pQCD now has a host of machinery dedicated to precision physics. Tools 
such as the Operator Product Expansion (OPE) are used to calculate many 
observables upto high orders in perturbation theory. Factorization theorems 
have been derived to separate out soft and hard dynamics in a systematic 
way; the former arc parameterized as non-pcrturbative parton distribution 
functions while the latter can be computed order by order in perturbation 
theory. The factorization of hard and soft scales is manifest in a renormal- 
isation group treatment that requires that physics be independent of this 
scale. An important result is the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) [8] equation for the evolution of quark and gluon distribution func- 
tions as a function of which made earlier work on the OPE accessible to 
straightforward experimental analysis. A consequence of DGLAP evolution 
is that, as shown in fig. 3, the parton distributions increase rapidly at small- 
x as increases. This means that as one probes finer and finer transverse 
resolution scales, more and more partons (which share the total hadron mo- 
mentum) can be resolved within the hadron. However, albeit the number of 
partons increases with increasing , the phase space density-thc number of 
partons/area/Q^-decreases and the proton becomes more and more dilute. 

1.2 The Regge-Gribov limit in QCD and Saturation 

The DGLAP framework in pQCD has worked very well and explains many 

features of the HERA data very well. However, when pushed to lower Q^, one 
begins to note increasingly unpleasant results. As shown in the right plot of 
fig. 3, fits to data push the gluon distributions into negative territory. This 
also seems to afflict Fl, which much be a positive definite quantity as it is 
directly proportional to a cross-section. Another feature of HERA data that 




the conventional pQCD approach needs further parameters to explain is the 
high fraction of diffractive events and the fact that these diffractive events have 
the same dependence as the total cross-section on the energy of the photon- 
hadron system. Within pQCD itself, it is expected that DGLAP should be 
supplemented by power corrections in called "higher twist" effects, which 
become increasingly important as is decreased and as one goes to smaller 
values of x. The higher twist formalism is however very cumbersome and not 
under theoretical control. It is therefore useful to take stock of the physics 
underlying the dynamics of partons in this kinematic region and approach 
the problem anew. 

The physics of the small x regime is best understood in a very different 
asymptotic limit from the better known (and understood) Bjorken limit and 
goes by the name of Regge-Gribov limit in QCD. This limit corresponds to 
a^Bj and s ^ oo with ^qcd) ~ fixed. The Regge-Gribov limit 

corresponds to the regime of strong color fields in QCD. It is responsible for the 
bulk of multi-particle production in QCD. At sufficiently large ^ ^qcD' 
the physics of this regime, albeit non-perturbative, should be accessible in 
weak coupling. 

In the Regge-Gribov limit, the leading logarithms in x, Q!5ln(l/x) dom- 
inate over the leading logs ols\vl{Q^ j A^fj^^-Xhe converse is true of course 
in the Bjorken limit. The renormalization group equation that resum these 
leading logs of x is called the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equa- 
tion [9] . The BFKL equation leads to a rapid power law growth in the gluon 
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distribution with decreasing x. A schematic diagram of the parton content of 
the hadron in the two asymptotic hmits is shown in fig. 4. The physics issues 




at small-x have a simple intuition. Both theory and experiments show that the 
gluon and the sea quark densities grow very rapidly in low-x. As parton dis- 
tributions grow, for fixed Q^, the occupation number of gluons in the hadron 
wavefunction becomes increasingly large and it is no longer possible to neglect 
their mutual interactions. In terms of Fcynman diagrams, these interactions 
involve gluon recombination and screening which deplete the gluon density 
relative to the bremsstrahlung diagrams that are responsible for the rapid 
growth of the parton distributions. The net effect of the competition between 
these two effects is that the occupation number of partons "saturates" at the 
maximum possible value oil/ as in QCD. The dynamics of QCD in this regime 
is fully non-linear corresponding to the strongest electric and magnetic fields 
in nature. Each mode in the wavefunction, because it has a different occupa- 
tion number, will saturate at a different value of x-the momentum scale Qs 
at which it saturates is called the saturation scale [10]. In terms of the gluon 
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distribution, the maximal phase space density is reached when 



1 xG{x,Ql) 



1 



(7) 



2(7V2 _ 1) ^^2Q2 



<^s{Ql) ' 



where Tri?^ is the hadron area in the impact parameter space which is only 
well defined if the wavelength of the probe is small compared to R, which we 
will assume throughout in the future discussions. As suggested by eq. 7, the 
growth of the gluon distribution functions with decreasing x implies a growth 
of the saturation scale as well. The saturation scale is the typical momen- 
tum of gluons in the high energy hadron wavefunction, and when x is small, 
Ql >• /Iqcd is a semi- hard scale in QCD. Thus as{Ql) << 1, suggesting that 
weak coupling techniques can be applied to study the saturation regime in 
the Regge-Gribov limit of QCD. The use of weak coupling techniques does 
not necessarily mean that the physics is perturbative; in this case, the many- 
particle interactions are responsible for the non-perturbativc nature of satura- 
tion, because, although individual interactions are small, the large number of 
gluons amplifies the effect necessitating a resummation of high parton density 
effects. 

To summarize, saturation is a natural consequence of QCD in the Regge- 
Gribov limit of the theory because gluon occupation numbers cannot be arbi- 
trarily large in the theory. Because saturation is achieved for different modes 
at different values of x, this also naturally suggests a dynamical scale that 
characterizes the onset of the non-linear saturation dynamics in QCD. If this 
scale is large compared to the QCD confining scale ~ ^qcd, the physics of 
saturation can be described using a weak coupling formalism. Because occupa- 
tion numbers are large in this regime, the appropriate degrees of freedom here 
arc classical fields. One is therefore, on very general grounds, led to postulate 
the existence of a weak coupling effective theory that captures the physics of 
this intensely non-linear regime of QCD. This is the Color Glass Condensate 
(CGC), and as we shall discuss, its name captures the properties of the matter 
in the hadron wavefunctions at high energies. 

1.3 The Color Glass Condensate 

One way to approach the small x problem is to appreciate that there is a 
formal Born-Oppenheimer separation between large x and small x modes [11] 
for quantum field theories on the light cone. These arc respectively the slow 
and fast modes in the effective theory. Thus on the time scale of the "wee" 
parton small x fields, the large x partons can be viewed as static charges. 
Since these are color charges, they cannot be integrated out of the theory but 
must be viewed as sources of color charge for the dynamical wee fields. With 
this dynamical principle in mind, one can write down an effective action for 
wee partons in QCD at high energies [12]. The generating functional of wee 
partons has the form 
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( 



(8) 



where the wee parton action has the form 



S[A,p] = ^ I 



d^xFl, + ^ [ d^x^dxSix-) X Tr {p{x^)U_^,^[A-]){9) 



In Eq. 8, p is a classical color charge density of the static sources and W[p] 
is a weight functional of sources (which sit at momenta fc+ > note, 
X = / P{^i^dron) ■ ^^"^ sources are coupled to the dynamical wee gluon fields 
(which in turn sit at < /l^) via the gauge invariant tcrm^ which is the 
second term on the RHS of Eq. 9. The first term in Eq. 9 is the QCD field 
strength tensor squared the wee gluons are treated in full generality in this 
effective theory, formulated in the light cone gauge j4+ = 0. The source j is an 
external source-derivatives taken with respect to this source (with this source 
then put to zero) generate correlation functions in the effective theory. 

The argument for why the sources are classical is subtle and follows from 
a coarse graining of the effective action to only include modes of interest. For 
large nuclei, or at small x, wee partons couple to a large number of sources. 
For a large nucleus, it can be shown explicitly that this source density is 
classical [14]. Further, it was conjectured that the weight functional for a large 
nucleus was a Gaussian in the source density (corresponding to tlw^ quadratic 
Casimir operator) [12, 15]. This was shown explicitly later to be correct albeit 
with a small but interesting correction, proportional to the cubic Casimir 
operator, that generates Odderon excitations in the effective theory [14]. For 
a large nucleus, the variance of the Gaussian distribution, the color charge 
squared per unit area p\, proportional to A^^^, is a large scalc-and is the 
only scale in the effective action'"'. Thus for >> ^qcd' '^s{Pa) << 1, and 
one can compute the properties of the theory in Eq. 8 in weak coupling. 

The Yang-Mills equations can be solved analytically to obtain the classi- 
cal field of the nucleus as a function of p: Aci.{p) [12]. From the generating 
functional in Eq. 8, one obtains for the two point correlator, 



Prom this expression, one can determine (for Gaussian sources) the occupation 
number (p = dN /-kE? /dk\dy of wee partons in the classical field of the nucleus. 
For A;_L >> Ql, one has the Weizsacker- Williams spectrum cj) Ql/k\, while 
for < Qs, one has a resummation to all orders in fc_L, which gives 4> ~ 
\n{Qs/k±). (The behavior at low k±_ can, more accurately, be represented 

* An alternative gauge invariant form of the coupling of sources and fields is ob- 
tained in Ref. [13]-it reproduces the BFKL equation more efficiently. 
PA is simply related to Qs'- Qs Q.6 pA- For a detailed discussion, see Ref. [16]. 




(10) 
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as ■^r{0, z) where F is the incomplete Gamma function and z = fc^/Q^.) A 
nice expression for the classical field of the nucleus containing these two limits 
is given in Ref. [17]. 

We are now in a position to discuss why a high energy hadron behaves 
like a Color Glass Condensate [5]. The "color" is obvious since the degrees of 
freedom, the partons, are colored. It is a glass because the stochastic sources 
(frozen on time scales much larger than the wee parton time scales) induce a 
stochastic (space-time dependent) coupling between the partons under quan- 
tum evolution (to be discussed in the next section)-this is analogous to a spin 
glass in condensed matter physics. Finally, the matter is a condensate since 
the wee partons have large occupation numbers (of order I /as) and have mo- 
menta peaked about Qs- Just as in actual condensates, the number of gluons, 
for a fixed configuration of sources, has a non-zero value and has a magnitude 
of order l/as- Gauge invariant obscrvablcs are computed by averaging over 
all possible configurations. As we will discuss, these properties are enhanced 
by quantum evolution in x. The classical field retains its structure-while the 
saturation scale grows: Qs{x') > Qs{x) for x' < x. 

The problem of small fluctuations about the effective action in Eq. 9 were 
first addressed in ref. [18] and it was noted that these gave large corrections of 
order 0:5 ln(l/.T) to the classical action; this implies that the Gaussian weight 
functional is fragile under quantum evolution of the sources. A Wilsonian 
renormalization group (RG) approach was later developed that systematically 
treated these corrections [19]. The basic recipe is as follows. Begin with the 
generating functional in eq. 8 at some with an initial source distribution 
W[p\. Perform small fluctuations about the classical saddle point of the effec- 
tive action, integrating out momentum modes in the region 
ensuring that A'^ is such that as ln(vl+/yl'^) << 1. The action reproduces 
itself at the new scale A'^ , albeit with a charge density p' = p + 5p, where Sp 
has support only in the window A'^ < fc+ < and [/?] — > Wjii+[p']. 
The change of the weight functional W[p] with x is described by the JIMWLK- 
non-linear RG equation [19] which we shall not write explicitly here-it can be 
found for instance in ref. [5, 6]. The JIMWLK equation has been re derived 
subsequently by several authors. We will discuss briefly in the next lecture, 
one such derivation by one of us and collaborators, in the context of nucleus- 
nucleus collisions. 

The JIMWLK equations form an infinite hierarchy (analogous to the 
BBGKY hierarchy in statistical mechanics) of ordinary differential equations 

for the gluon correlators < A1A2 ■ ■ ■ An >y, where Y = ln(l/.x) is the rapidity. 
The expectation value of such an operator O is defined to be 

{0)y= j[da]0[a\WY[a], (11) 

where a = p. The corresponding JIMWLK equation for this operator is 
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d{0[a])Y ,1 / 5 , u 



Here x here is a non-local object expressed in terms of path ordered (in ra- 
pidity) Wilson lines of a [5]. This equation is analogous to a (generalized) 
functional Fokker-Planck equation, where Y is the "time" and x is a general- 
ized diffusion coefficient. This equation illustrates the stochastic properties of 
operators in the space of gauge fields at high energies. For the gluon density, 
which is proportional to a two-point function < a"(x_L)a^(y_L) >, one recovers 
the BFKL equation in the limit of low parton densities. For a first attempt to 
solve the JIMWLK equation numerically, sec rcf. [20]. 

For large Nc and large A {ogA^/'^ » 1), the expectation value of the 
product of traces of Wilson lines factorizes into the product of the expectation 
values of the traces: 

< tv{v,v})Tt{v,v^) >^< tv(\4f;) > < TV(v;v;t) > , (13) 

where Vx = 'Pexp(J dz~a°'{z~ ,x_[_)T°'). Here V denotes path ordering in 
x~ and is the SU(3) generator in the adjoint representation. In the dipole 
picture, the cross-section for a dipole scattering off a target P can be expressed 
in terms of these 2-point dipole operators as [36] 



(Tqqp{x,rA_) =2 j (fbAfY{x,r±,b), 



(14) 



where My, the imaginary part of the forward scattering amplitude, is defined 
to be A/y = 1^77-(Tr(T4yjJ))y. Note that the size of the dipole, r± = x±—y^ 
and b = {x±_ + y_L)/2. The JIMWLK equation for the two point Wilson 
correlator is identical in this large A, large A^c mean field limit to an equation 
derived independently by Balitsky and Kovchegov-the BK equation [21], which 
has the operator form 

^ = as /Cbfkl (^{My-J^] . (15) 

Here /Cbfkl is the well known BFKL kernel. When M « I, the quadratic 

term is small and one has BFKL growth of the number of dipoles; when A/" is 
close to unity, the growth saturates. The approach to unity can be computed 
analytically [22]. The BK equation is the simplest equation including both 
the Brcmsstrahlung responsible for the rapid growth of amplitudes at small 
x as well as the repulsive many body effects that lead to a saturation of this 
growth. 

In this framework, a saturation condition, say M = 1/2, determines the 
saturation scale. One obtains = Qq exp(AF), where A = cas with c « 4.8. 
As we shall discuss further in the next subsection, the (arbitrary) choice of 
this saturation condition affects the overall normalization of this scale but 
does not affect the power A. In fixed coupling, the power A is large and 
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there are large pre-asymptotic corrections to this relation-which die off only 

slowly as a function of Y. BFKL running coupling effects change the be- 
havior of the saturation scale completely-one goes smoothly at large Y to 
Q2 = Qlcxp{^2bociY + Yo)) where 60 is the coefRcient of the one-loop QCD 
/3-function. An impressive computation of Q,, is the work of Triantafyllopou- 
los, who obtained Qs by solving NLO-resummed BFKL in the presence of an 
absorptive boundary (which approximates the CGC) [23]. The pre-asymptotic 
effects are much smaller in this case and the coefficient A « 0.25 is very close 
to the value extracted from saturation model fits to the HERA data. There is 
currently an intense, on-going effort to directly compute the NLO corrections 
to the leading order kernels of the BK equation [27]. No analytical solution 
of the leading order BK equation exists in the entire kinematic region but 
there have been several numerical studies at both fixed and running cou- 
pling [24, 25, 26]. These studies suggest that the solutions have a soliton like 
structure and that the saturation scale has the behavior discussed here. 

The soliton like structure of the numerical solutions is not accidental, as 
was discovered by Munier and Peschanski [28]. They noticed that the BK- 
equation, in a diffusion approximation, bore a formal analogy to the FKPP 
equation describing the propagation of unstable non-linear wavefronts in sta- 
tistical mechanics [29]. In addition, the full BK-cquation lies in the universality 
class of the FKPP equation. This enables one to extract the universal proper- 
ties such as the leading pre-asymptotic terms in the expression for the satu- 
ration scale. It was realized [30] that a stochastic generalization of the FKPP 
equation-the sFKPP equation-could provide insights into impact parameter 
dependent fluctuations [31] in high energy QCD beyond the BK-equation. We 
shall not discuss further efforts in that direction here except to note that there 
are many open ends (and opportunities) here which are still not satisfactorily 
resolved and require new ideas. 

1.4 Color Dipole CGC models in DIS and hadronic scattering 

In the previous sub-section, we discussed the formalism of the CGC and the 
JIMWLK/BK RG equations. We will here discuss phenomenological appli- 
cations of these ideas to DIS and (more briefly) to hadronic collisions. The 
next lecture will discuss more at length the application of these ideas to heavy 
ion collisions. A strong hint that semi-hard scales may play a role in small 
X dynamics came from "geometrical scaling" of the HERA data [4]. The in- 
clusive virtual photon-|-proton cross-section for x < 0.01 and all available 
scales ^ as a function of r = Q^/Ql, where Ql{x) = exp(AF) GeV^. Here 
Y = ln{xo/x) is the rapidity; xq = 3 ■ 10~^ and A = 0.288 are parameters 
fit to the data [4, 32]. Geometrical scaling of the e-|-p data is shown in fig. 6, 
which demonstrates that the inclusive diffractive, vector meson and DVCS 



® The E665 data are a notable exception. 
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Fig. 5. Interaction of the color dipole with the nuclei 

cross-sections at HERA, with a slight modification^in the definition of r, also 
appear to show geometrical scaling [32]. A recent "quality factor" statisti- 
cal analysis [33] indicates that this scaling is robust; it is however unable to 
distinguish between the above fixed coupling energy dependence of Qs and 
the running coupling Qs{x) oc exp(\/y) dependence of the saturation scale. 
Geometrical scaling is only asymptotic in both fixed and running coupling 




Fig. 6. Geometrical scaling fully inclusive, diffractive and exclusive vector meson 
cross-sections. From [32]. 

evolution equations. However, recent analyses [34, 35] suggests that the onset 
of the scaling asymptotics may be precocious. Geometrical scaling alone is not 
sufficient evidence of saturation effects and it is important to look at the data 
in greater detail in saturation/CGC models. 



^ td,vm = iQ^ + M'^)/Q'1, where M denotes the mass of the diffractive/ vector 
meson final state. 
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All saturation models [36] express the inclusive virtual photon+proton 
cross-section as 



Here 



(Tf = d rj_ dz 



^2.j,{r_L,z,Q) 



2 



^^ZtI^-L-Q.^) • <^qqPi'r±,XBi,b±) (16) 



represents the probability for a virtual photon to pro- 
duce a quark - anti-quark pair of size r = \r±\ and (T(r^, .tbj, denotes the 
dipole cross section for this pair to scatter off the target at an impact param- 
eter b±. The former is well known from QED, while the latter represents the 
dynamics of QCD scattering at small x see eq. 14. A simple saturation model 
(known as the GBW model [37]) of the dipole cross section, parametrized 
as aqqp = 2(1 — e~'"^'5^'f (^V4) where Qs,p{x) = (xq/x)"^ GeV^, gives a good 
qualitative fit to the HERA inclusive cross section data for xq = 3 ■ 10^^ and 
A = 0.288. Though this model captures the qualitative features of saturation, 
it does not contain the bremsstrahlung limit of perturbative QCD (pQCD) 
that applies to small dipoles of size r <^ ^/Qs.p{x)- 

In the classical effective theory of the CGC, one can derive, to leading 
logarithmic accuracy, the dipole cross section [38] containing the right small 
r limit. This dipole cross section can be represented as [39] 

(Tqqp{r±,XBi,b±) = 2 [1 - exp {-r''F{xBi,r±)Tp{b±))] , (17) 

where Tp{b±) is the impact parameter profile function in the proton, normal- 
ized as J d^bTp{b±) = 1 and F is proportional to the gluon distribution [40]: 

FixBi,rl) = n^as (/x^ + 4/ri) xb^Q (a;Bj,Mo + /(2A^c) , (18) 

evolved from the initial scale Ho by the DGLAP equations. The dipole cross 
section in eqn. 17 was implemented in the impact parameter saturation model 
(IPsat) [39] where the parameters are fit to reproduce the HERA data on 
the inclusive structure function F2. Here p is defined^ as the solution of 

(^qqp{rx,XBi,b±) = I / Q sA^Bj, bx) =2(1-6"^/'*). 

The IPsat dipole cross section in eqn. 17 is valid when leading logarithms 
in X in pQCD are not dominant over leading logs in Q^. At very small x, 
where logs in x dominate, quantum evolution in the CGC describes both the 
BFKL limit of linear small x evolution as well as nonlinear JIMWLK/BK 
evolution at high parton densities [19, 21]. These asymptotics are combined 
with a more realistic 6-dependence in the b-CGC model [41, 42]. Both the 
IPsat model and the b-CGC model provide excellent fits to HERA data for 
X < 0.01 [42, 43]. The saturation scale extracted from the fit in the IPsat 
model is shown in Fig. 7. The important point to note is that the energy 
dependence of the extracted Qs,p is significantly stronger than those predicted 
in non-perturbative models [44]. 



* This choice of is ccjuivalcnt to the saturation scale in the GBW model for the case 
of a Gaussian dipole cross section. 
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Fig. 7. The saturation scale vs Ijx in the IPsat and b-CGC models [42]. 



The strong field dynamics of small x partons is universal and should be 
manifest in large nuclei at lower energies than in the proton. In Fig. 8 (left), 
we show the well known shadowing of in the fixed target e+A E665 and 
NMC experiments. Expressed in terms of r = /Ql (Fig. 8 (right)), the data 
show geometrical scaling [45] . A careful study of nuclear DIS in the IPsat CGC 



2 F "/AF, 



0.4^ 
10' 
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10"' II 



Fig. 8. Left: Shadowing of F2 from the NMC and E665 fixed target experiments. 
Right: The data scaled as a function of r = /Qi,A [45]. 



framework was performed in Ref. [39, 46]. The average differential dipole cross 
section is well approximated by {■j§^)N ~ 



TA{b±). 



where 



Ta(^j_) is the well known Woods Saxon distribution. The average is defined 
as {0)n = JYli^i(Pb^ iTA{b±j)0{{bj^^i}). Here cTp is determined from the 
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IPsat fits to the e+p data; no additional parameters are introduced for eA 
collisions. In Fig. 9 (left), the model is compared to NMC data on Carbon 
and Calcium nuclei-the agreement is quite good. In Fig. 9 (right), we show 
the extracted saturation scale in nuclei for both central and median impact 
parameters. To a good approximation, the saturation scale in nuclei scales as 
QlAix,bmed.) ~ Qlp{x,b^ed.) ■ (A/x)^/^. The factor of 2001/^ « 6 gives a 
huge "oomph" in the parton density of a nucleus relative to that of a proton 
at the same x. Indeed, one would require a center of mass energy ~ 14 times 
larger in the proton case. At extremely high energies, this statement must be 
qualified to account for the effects of QCD evolution [47]. 



Fig. 9. Left: Comparison of the IPsat model (with no adjustable parameters) to 
the NMC data. Right: The A and x dependence of the saturation scale in the IPsat 
model [46]. 

We now turn to a discussion of CGC effects in hadronic collisions. A sys- 
tematic treatment of the scattering of two strong color sources (such as two 
high energy nuclei) is discussed in the next lecture. To leading order, the 
problem reduces to the solution of the classical Yang-Mills (CYM) equations 
averaged over color sources for each nucleus [48, 49]; the variance of this dis- 
tribution of sources is proportional to ^. Besides the nuclear radius, it is 
the only scale in the problem, and the ^ ~ ■ {A/x)^'^ expression for the 
saturation scale was used in CGC models of nuclear collisions to successfully 
predict the multiplicity [49] and the centrality dependence of the multiplic- 
ity [73] dependence in gold-|-gold collisions at RHIC. The universality of the 
saturation scale also has a bearing on the hydrodynamics of the Quark Gluon 
Plasma (QGP); the universal form leads to a lower eccentricity [50] (and 
therefore lower viscosity) than a non-universal form that generates a larger 
eccentricity [51] (leaving room for a larger viscosity) of the QGP^ 

For asymmetric (off-central rapidity) nuclear collisions, or proton/deuteron 
-|- heavy nucleus collisions, fc_L-factorization can be derived systematically for 
gluon production, at leading order, in the CGC framework [52]. Limiting frag- 

^ We shall present yet another take on this issue in lecture II. 




oc 
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mentation [53] and deviations thereof, are described by solutions of the BK- 
equation. Predictions for the multiphcity distribution in A+A colhsions at the 
LHC [54] for both Golec-Biernat-Wusthoff (GBW) and classical CGC (MV) 
dipole initial conditions give a charged particle multiplicity of 1000-1400 in 
central lead+lead collisions at the LHC^^. The results are shown in Fig. 10. 




n Ibeam 



Fig. 10. Prediction for limiting fragmentation and deviations away from it at LHC 
energies. The bands denote the range in the predictions for GBW and MV models. 
From [54] . 

In deuteron + gold collisions at RHIC, data on the inclusive hadron spec- 
tra^^ can be directly compared to model predictions [55] . The result is shown 
in Fig. 11. For a comprehensive review of applications of CGC picture to 
RHIC phenomenology, we refer the reader to Ref. [57]. There are a couple 
of caveats to this picture. Firstly, k± factorization is very fragile. It does not 
hold for quark production even at leading order in the parton density [58], 
albeit it may be a good approximation for large masses and transverse mo- 
menta [59]. For gluon production, it does not hold beyond leading order in 
the parton density [60, 49]. Secondly, a combined comprehensive analysis of 
HERA and RHIC data is still lacking though there have been first attempts 
in this direction [61]. 

1.5 The future of small x physics at hadron colUders and DIS 

The LHC is the ultimate small x machine in terms of reach in x for large . 
A plot from Ref. [62] illustrating this reach is shown in Fig. 12 (left). For a 

^° The MV initial condition has the same form as the IPsat dipole cross-section 
discussed earlier. 

See Ref. [74] for other model predictions. 

The same analysis also gives good agreement for the forward p+p spectrum at 
RHIC [56]. 
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Fig. 11. The inclusive k± distributions in deuteron+gold collisions compared to 
theory curves for different rapidities. From [55]. 



recent review of the small x opportunities at the LHC, see Ref. [63]. The LHC 
will provide further, more extensive tests of the hints for the CGC seen at 
RHIC. The universality of parton distributions is often taken for granted but 
factorization theorems proving this universality have been proven only for a 
limited number of inclusive final states. However, as we have discussed, small 
X is the domain of rich multi-parton correlations. These are more sensitive to 
more exclusive final states for which universality is not proven [64]. Therefore, 
while the LHC will have unprecedented reach in x, precision studies of high 
energy QCD and clean theoretical interpretations of these motivate future 
DIS projects. Two such projects are the EIC project in the United States [65] 
and the LHeC project in Europe [66]. As we discussed previously, strong color 




Fig. 12. Left: Kinematic x-Q^ reach of different final states at the LHC compared 
to other experiments with nuclei . From [62]. Right: The saturation scale in the 
proton, calcium and gold in the liinematic acceptance of the EIC. 
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fields may be more easily accessible in DIS off nuclei relative to the proton due 
to the "oomph" factor. In Fig. 12 (right), we show the saturation scale Qs,a{x) 
overlaid on the x-Q^ kinematic domain spanned by the ETC. As suggested by 
the figure, the EIC (and clearly the higher energy LHeC) will cleanly probe 
the cross-over regime from linear to non-linear dynamics in QCD. For further 
discussion of the physics of an Electron Ion collider, see Ref. [67, 68]. 

Heavy Ion Collisions 

In the first lecture, we discussed the physics of high parton densities in QCD. 
We motivated a description of the high energy structure of hadrons and nuclei 
as a Color Glass Condensate. We discussed the phenomenological application 
of CGC based models to describe data on DIS and hadronic collisions. In this 
lecture, we will focus our attention on heavy ion collisions and try and un- 
derstand ab initio the properties of the QCD matter that is formed when two 
high energy nuclear wavefunctions-sheets of colored glass-collide. A schematic 
space-time picture of the evolution of a heavy ion collision is shown in fig. 13. 

We shall not attempt to describe all features of a heavy ion collision but 
merely focus on the very early initial dynamics for its intrinsic interest but 
also because it is important to understand to determine whether the matter 
thermalizes and can be subsequently described by hydrodynamics. The QCD 
matter that is formed at very early times is a coherent classical field, which 
expands, decays into nearly on shell partons and may eventually thermalize 
to form a Quark Gluon Plasma (QGP). Because it is formed by melting the 
frozen CGC degrees of freedom, and because it is the non-equilibrium matter 
preceding the QGP, this matter is called the Glasma [69]. 



t 




hydrodynamics 



classical EOMs 



kinetic theory 



Fig. 13. The various stages of a heavy ion collision. 
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2 From CGC to Glasma 



We will first discuss the classical picture of the Glasma that emerges from 
solutions of classical Yang-Mills equations. We will then discuss the role of 
quantum corrections in the Glasma. 

2.1 The Classical Solution 

Let us begin by setting up the kinematics involved for the problem of two ultra- 
relativistic nuclei approaching each other. They are highly Lorcntz contracted 
and can be considered to be sitting at = for t < in the light cone 
coordinates. They collide at a; = < = 0. At i > it is convenient to introduce 
the proper time r = \Jt'^ — which is invariant under Lorentz boosts and 
the space time rapidity 



For free streaming particles z = vt ~ pzt/ E and the space-time rapidity equals 
the momentum space rapidity y, namely, rj = y. 



Fig. 14. The space-time distribution of gauge fields before and after the colhsion. 

Let us make some some approximations to get to the solutions without 
sacrificing any essential physics. The collisions are considered to be at very 
high energy so that the nuclei are static light cone currents that are delta 
functions in respectively for nuclei whose large momentum component is 
given by P^. A consequence is that the solutions of the classical Yang-Mills 
equations, for this source distribution, are boost invariant. This leads to a 
considerable simplification since the equations now only depend on the two 




(19) 




pure gauge 2 
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Fig. 15. The color electric and the magnetic fields lying on the sheets of the charge 
distribution. They are perpendicular to each other and to the beam direction. 

transverse directions and the proper time. The initial conditions for the evo- 
lution of the classical gauge fields in the forward light cone is illustrated in fig. 
14. The fields of the nuclei before the collision are the Lienard-Wiechert po- 
tentials associated with the respective color charge densities. Both the charge 
densities and the fields only exist on the sheets and for each source of charge 
the electric field E and the magnetic field B are orthogonal to each other and 
to the beam direction. The situation before the collision is depicted in fig 15. 
However the vector potentials lie outside the sheets but possess a discontinuity 
across the sheet according to Gauss law. This is where the small-x gluons are 
located having a large longitudinal extent and coupling with a host of color 
sources. As has been argued before they are represented by the classical color 
fields, have very small lifetime and see the color sources to be static during 
their lifetime. As noted previously, the infinitesimal nature of the sheets is 
intended to simplify the problem. A finite size can be taken care of using RG 
arguments. Fig. 14 shows that the fields vanish in the backward light cone and 
are two dimensional pure gauge transforms of vacuum in both nuclei before 
the collision. 

In the forward light cone, a pure gauge solution of the Yang-Mills (YM) 
equations of motion cannot be found-the sum of two pure gauges in QGD 
is not a pure gauge. Solving the YM equations near the light cone [49, 72, 
70], with the pure gauge initial conditions from the two nuclei before the 
collision, one finds that the transverse color E and B fields vanish as r ^ 
but the longitudinal fields are non- vanishing. The results of a numerical 
simulation [70] are shown in fig. 16. We note that the non-zero E and B fields 
imply that the initial conditions [71] may have a large density of topological 
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charge F'^^F^u- The LO picture of the Glasma that emerges is one of color 
flux tubes carrying non-trivial topological charge, localized in the transverse 
plane and of transverse size l/Q^, stretching between the valence color degrees 
of freedom. As we shall soon see, this picture provides a plausible explanation 
of the near side ridge in heavy ion collisions. 

For the inclusive gluon distribution, the leading order (LO) classical con- 
tribution is of order 0{-^) but all orders in gpi.2, where pi_2 are the local 
color charge densities of the two nuclei. It can be expressed as 

xY^e'iipy^^p) {A,{x)A,{y)) . (20) 

A 

The gauge fields on the right hand side are computed numerically for proper 
times r > [72, 70] by solving the classical Yang-Mills (CYM) equations in the 
presence of the light cone current J^^'°^ = 5^^ 5{x~)p1{x±) + 5{x^)p2{x±^. 
The expectation value (• • •) for any inclusive operator O denotes^^ 

{O) = [ [Dpi] [Dp,] _y [pi] [h] O [pi ,P2] ■ (21) 



Here p is the local color charge density in Lorentz gauge, which is related by a 
simple gauge transformation to the corresponding color charge density p in light 
cone gauge. We note that W[p] is gauge invariant. The light cone gauge classical 
fields are expressed most simply in terms of color charge densities in Lorentz 
gauge. 



24 Debasish Banerjee, Jajati K. Nayak, and Raju Venugopalan 

We will justify the validity of this formula shortly. For central rapidity RHIC 

heavy ion collisions, as discussed in the previous lecture, evolution effects arc 
not important, and VF[/o] is a Gaussian functional in p with the variance 
Recall that Qs ~ 0.6 /xa- So, for these Gaussian distributions, performing the 
average in eq. 21 over the solutions to the Yang-Mills equations in eq. 20, 
one can compute the number (and energy distributions) of produced gluons 
in terms of Qs- For the energy density, one obtains £ 20 - 40 GeV/fm^ for 
the values of Qs we mentioned — obtained by extrapolating from fits to the 
HERA and fixed target e+A data. 

This LO formalism was applied to successfully predict the RHIC multi- 
plicity at y ^ [72, 70] as well as the rapidity and centrality distribution 
of the multiplicities [73]. At LO, the initial transverse energy is Et ~ Qs, 
which is about 3 times larger than the final measured Et, while (assuming 
parton hadron duality) iVcGC ^ -^had.- The two conditions arc consistent if 
one assumes nearly isentropic flow which reduces Et due to PdV work while 
conserving entropy. This assumption has been implemented directly in ideal 
hydrodynamic simulations [75]. 

As discussed previously, CGC based models give values for the initial ec- 
centricity e that are large than those in Glauber model because the energy 
and number density locally is sensitive to the lower of the two saturation 
scales (or local participant density) in the former and the average of the two 
in the latter. Naively, CGC initial conditions would have more flow then and 
have more room for dissipative effects relative to Glauber. This conclusion is 
turned on its head in a simple parameterization of incompletely thermalized 
flow [76]: V2/€ = jy+RJk^^ where K = ■^<^^Cs is the Knudsen number, a 
the cross-section, Cg the sound speed and S± the transverse overlap area. Kq 
is a number of order unity. If thermalization were complete, K ^ Q and one 
approaches the hydro bound. Computing e with different initial conditions, 
and plotting the l.h.s ratio versus one has a two parameter flt to a 

and Cs- The greater CGC eccentricity forces W2/e to be lower for more central 
collisions thereby leading to lower Cg; quicker saturation of V2/e forces larger 
a and therefore lower rj in the CGC relative to Glauber [77]. While is con- 
ceivable however that this result may not prove robust against more detailed 
modeling, it is clear that the results are very sensitive to the initial conditions. 

How much flow is generated in the Glasma before thermalization? The 
primordial Glasma has occupation numbers / ~ ^ and can be described as 
a classical field. As the Glasma expands, higher momentum modes increasingly 
become particle like and eventually the modes have occupation numbers / < 1, 
which may be described by a thermal spectrum. A first computation of elliptic 
flow of the Glasma found only about half the observed elliptic flow [78] albeit 
the computation did not properly treat the interaction between hard and 
soft modes in the Glasma. Formulating a kinetic theory that describes this 
evolution is a challenging problem in heavy ion coUisions-for a preliminary 
discussion, see [79]. 
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Fig. 17. Left: Chern-Sirnons mean squared charge as a function of proper time 
(in units of the saturation scale), generated in the LO boost-invariant 2-l-l-D 
Glasma [71]. Quantum fluctuations allow sphaleron transitions which may give rise 
to significantly greater values of \/Tp^- Right; The same quantum fluctuations, al- 
beit suppressed by as, grow as as exp(v^Qsr) [80]. These unstable fluctuations, 
resummed to all orders, likely give rise to more isotropic initial distributions in the 
Glasma [87]. 

The LO field configurations are unstable and lead to very anisotropic mo- 
mentum distributions at later times. Such distributions can trigger an insta- 
bility analogous to the Weibel instability in QED plasmas [81]. It is observed 
in 3-1- 1-D numerical solutions of GYM equations [80] that small rapidity de- 
pendent quantum fluctuations in the initial conditions generate transverse E 
and B fields that grow rapidly as exp(-\/Q5r). They are the same size as the 
rapidly diluting longitudinal E and B fields on time scales of order ^ In^(^). 
The transverse E and B fields may cause large angle deflections of colored 
particles leading to -p^ broadening and energy loss of jets-recent numerical 
simulations by the Frankfurt group appear to confirm this picture [82] . These 
interactions of colored high momentum particle like modes with the soft co- 
herent classical field modes may also generate a small "anomalous viscosity" 
whose effects on transport in the Glasma may mask a larger kinetic viscos- 
ity [83]. The same underlying physics may cause "turbulent isotropization" 
by rapidly transferring momentum from soft "infrared" longitudinal modes 
to ultraviolet modes [84]. Finally, albeit the LO result demonstrated that 
one could have non-trivial Ghern-Simons charge in heavy ion collisions, the 
boost invariance of GYM equations disallows sphaleron transitions that per- 
mit large changes in the Ghern-Simons number [71]. With rapidity dependent 
quantum fiuctuations, sphaleron transitions can go. These may have impor- 
tant consequences-in particular P and CP odd metastable transitions that 
cause a novel "Ghiral Magnetic Effect" [85] in heavy ion collisions. Numerical 
GYM results for Ghern-Simons charge and (square root) exponential growth 
of instabilities are shown in fig. 17. 
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2.2 QCD Factorization and the Glasma instability 

The discussion at the end of the last subsection strongly suggests that next- 
to-leading order (NLO) quantum fluctuations in the Glasma, while paramet- 
rically suppressed, may alter our understanding of heavy ion collisions in a 
fundamental way. To cosmologists, this will not come as a surprise-quantum 
fluctuations play a central role there as well. In recent papers, it was shown for 
a scalar theory that moments of the multiplicity distribution at NLO in A+A 
collisions could be computed as an initial value problem with retarded bound- 
ary conditions [86]; this framework has now been extended to QCD [87, 88]. 
In QCD, the problem is subtle because some quantum fluctuations occur 
in the nuclear wavefunctions and are responsible for how the wavefunctions 
evolve with energy; others contribute to particle production at NLO. Fig. 18 
illustrates particle production in fields theories with strong sources and the 
non-factorizable quantum fluctuations that are suppressed in the leading log 
framework. 

In writing eq. 21, its scope of validity was not specified. A factorization 
theorem organizing quantum fluctuations shows that all order leading loga- 
rithmic contributions to an inclusive gluon operator O in the Glasma can be 
expressed as [87, 88] 



(^^>L.o, = / \DmDp2\ M^y._-y[pi] W%_+y[p2] 0,o [?i,P2] , (22) 



where 0^0 the same operator evaluated at LO by solving classical Yang- 
Mills equations and VFYi,eamT>'[/'i>2] a-rc the weight functionals that obey the 
JIMWLK Hamiltonians describing the rapidity evolution of the projectile and 
target wavefunctions respectively. This theorem is valid^^ if the rapidity in- 
terval corresponding to the production of the final state, Ziy < The 
are analogous to the parton distribution functions in collincar factorization; 
determined non-perturbatively at some initial scale Ioj their evolution with 
Y is given by the JIMWLK Hamiltonian. 

This factorization theorem is a necessary first step before a full NLO com- 
putation of gluon production in the Glasma. Eq. (22) includes only the NLO 
terms that are enhanced by a large logarithm of l/xi,2, while the complete 
NLO calculation would also include non enhanced terms. These would be of 
the same order in as as the production of quark-antiquark pairs [89] from 
the classical field. To be really useful, this complete NLO calculation likely 
has to be promoted to a Next-to-Leading Log (NLL) result by resumming all 
the terms in qs(q:s lii(l/a;i,2))". Now that evolution equations in the dense 
regime are becoming available at NLO, work in this direction is a promising 
prospect. In addressing the role of instabilities at NLO, note that small field 
fluctuations fall into three classes: i) Zero modes (p,, = 0) that generate the 



Recent work suggests that this factorization theorem is valid even when the rar- 
pidity restriction is relaxed. 
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Fig. 18. Left: Illustration of particle production in a field theory with strong time 
dependent sources. Right: Quantum fluctuations, where the nuclei talk to each other 
before the collision, are suppressed; this is responsible for high energy factorization 
of inclusive gluon operators O in A+A collisions [87, 88]. 

leading logs resummed in eq. 22; the coefficients of the leading logs do not 
depend on . ii) Zero modes that do not contribute at leading log because 
they are less singular than the leading log contributions. These become rele- 
vant in resumming NLL corrections to the factorization result. Because they 
are zero modes, they do not trigger plasma instabilities, iii) Non zero modes 
{jJr, 7^ 0) that do not contribute large logarithms of l/a;i_2, but grow exponen- 
tially as exp(-\/Q5r). While these boost non-invariant terms are suppressed by 
as, they are enhanced by exponentials of the proper time after the collision. 
These leading temporal divergences can be resummed and the expression for 
inclusive gluon operators in the Glasma revised to read 

(0>LL.g+Li„., = / \D'kW>n\ H'V..„-rl?,l»'^..„,+vp2l 

x/[D„,.)]Z|aM10.„|.47+a,*+a| (23) 

where A^{x) = —-^pi{x±,x^) and yl^(a;) = —-gTP2{x±,x'^). The effect of 
the resummation of instabilities is therefore to add fluctuations to the initial 
conditions of the classical field, with a distribution that depends on the out- 
come of the resummation. This spectrum Z is the final incomplete step 
in determining all the leading singular contributions to particle production in 
the Glasma-see however Ref. [90] . The stress-energy tensor T^"^ can then be de- 
termined ab initio and matched smoothly to kinetic theory or hydrodynamics 
at late times. 

2.3 Two particle correlations in the Glasma and the Ridge 

Striking "ridge" events were revealed in studies of the near side spectrum of 
correlated pairs of hadrons at RHIC [91]. The spectrum of correlated pairs on 
the near side of the STAR detector extends across the entire detector accep- 
tance in pseudo-rapidity of order Arj ^ 2 units but is strongly coUimated for 
azimuthal angles Ac/). Preliminary analysis of measurements by the PHENIX 



28 Debasish Banerjee, Jajati K. Nayak, and Raju Venugopalan 



and PHOBOS collaborations corroborate the STAR results. In the latter case, 
the ridge is observed to span the PHOBOS acceptance in pseudo-rapidity of 
At] 4 units. 




Fie. 19. LeftiGlasma flux tubes of transverse size 7=7- < -r^ — . The field lines 

" QS ^QCD 

correspond to parallel E and B fields, which carry topological charge. Right: The 
leading two particle contribution. Superficially disconnected, they are connected by 
averaging over the large x color sources. Systematic power counting shows these 
graphs dominate over usual "pQCD" graphs at high energies. 

Causality dictates (in strong analogy to CMB superhorizon fluctuations) 
that long range rapidity correlations causing the ridge must have occurred 
at proper times r < Tfreeze out e~ 2 I^a"?'!? I ^ where and are the rapidi- 
ties of the correlated particles. If the ridge span in psuedo-rapidity is large, 
these correlations must have originated in the Glasma. The PHOBOS data 
suggest correlations at times as early as a fermi. As noted previously, parti- 
cles produced from Glasma flux tubes are boost invariant. See fig. 19. The 
correlated two gluon production in the Glasma flux tube is independent of 
rapidity [92]-thereby allowing long range correlations. Remarkably, the result 
shown in fig. 19 (right), that "classical" disconnected graphs give the leading 
contribution to two particle correlations, holds true even when one includes 
leading logarithmic contributions to all orders in perturbation theory [88]. 
When the separation in rapidity between pairs Ay 3> l/as; particle emission 
between the triggered pairs will modify this result. This rapidity dependence 
can be computed and tested at the LHC. 

Ours is the only dynamical model with this feature-for other models, see 
Ref.[93]. The particles produced in a flux tube are isotropic locally in the 
rest frame but are coUimated in azimuthal angle when boosted by transverse 
flow [94]. Combining our dynamical calculation of two particle correlations 
with a simple "blast wave" model of transverse flow, we obtained reasonable 
agreement with 200 GeV STAR data on the amplitude of the correlated two 
particle spectrum relative to the number of binary collisions per participant 
pair. A more sophisticated recent treatment of the flow of Glasma flux tubes 
shows excellent agreement of the amplitude of the ridge with centrality and 
energy, and likewise of the angular width of the ridge with centrality and en- 
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ergy [95]. Three particle correlations can further help distinguish the Glasma 
explanation for the ridge from other mechanisms. 
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